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The first dedicated search for ultra-high energy (UHE) tau neutrinos of astrophysical origin was 
performed using the IceCube detector in its 22-string configuration. The search also had sensitivity 
to UHE electron and muon neutrinos. After application of all selection criteria to approximately 
200 live-days of data, we expect a background of 0.60 ± 0.19 (stat.) t-o'gs i s Y s ) events and observe 
three events, which after inspection emerge as being compatible with background but are kept in 
the final sample. Therefore, we set an upper limit on neutrinos of all flavors from UHE astrophysical 
sources at 90% CL of E 2 &(v x ) < 16.2 x 10 -8 GeV cm" 2 sr" 1 s" 1 over an estimated primary neutrino 
energy range of 340 TeV to 200 PeV. 

PACS numbers: 95.85.Ry, 14.60.Lni, 95.30.Cq, 95.55.Vj, 14.60.Fg 
Keywords: IceCube, neutrino telescope, tau neutrinos, double bangs 



I. INTRODUCTION 

Proposed astrophysical sources of observed ultra-high 
energy (UHE) cosmic rays are expected to also produce 
ultra-high energy neutrinos, mainly via charged pion de- 
cay following interactions on ambient matter and radi- 
ation 0, [2] • Candidate neutrino sources include active 
galactic nuclei, gamma ray bursts and microquasars 0- 
Lg] . Neutrinos are expected to arrive at Earth with a flavor 
ratio of v c :v ll :v T = 1:1:1 in the standard neutrino oscilla- 



* Physics Department, South Dakota School of Mines and Tech- 
nology, Rapid City, SD 57701, USA 

t Los Alamos National Laboratory, Los Alamos, NM 87545, USA 

* also Sezione INFN, Dipartimento di Fisica, 1-70126, Bari, Italy 
§ Corresponding author: Seon-Hee Seo. seo@fysik.su.se 

1 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 



tion scenario [7|. Other neutrino production and propa- 
gation models predict different flux ratios at Earth [8 10] . 
If there are many astrophysical point sources of neutri- 
nos, but each one is too weak to be distinguished individ- 
ually from background, then a suitable detection strat- 
egy is to perform a cumulative search for "diffuse" flux 
of UHE neutrinos over the full available solid angle. 

In previous searches [TTI - fT3j , diffuse astrophysical UHE 
neutrinos were distinguished from atmospheric neutrinos 
by requiring the energy of candidate UHE neutrino events 
to exceed a certain threshold. In this work, we present 
techniques for identifying u T interactions and show the 
results of the first search for diffuse astrophysical UHE 
neutrinos that specifically selected events consistent with 
several v T interaction topologies. 

At E„ > 1 PeV, a search for UHE v T offers several 
advantages over UHE v e and f M searches. Partially and 
fully contained interactions of UHE v T in the detector can 
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FIG. 1. (Color online) A simulated double bang event with 
a primary v T energy of 47 PeV entering the IC22 detector at 
a zenith angle of 35°. The two bigger circles in gray color 
represent the vertexes of tau creation (upper left) and decay 
(lower right) which are connected by a tau track. Each smaller 
circle represents a photomultiplier tube (PMT) that observed 
light produced by the event. The arrival times are indicated 
by colors that are ordered following the colors of the rainbow, 
with red corresponding to earlier times and violet to later 
times. 



produce very distinctive signatures owing to the macro- 
scopic r decay length. Each such signature should ap- 
pear in proportion to the well-measured r branching ra- 
tios [HI, providing a useful cross-check on the positive 
identification of multiple v T events. As shown below, 
the chief sources of possible background events are un- 
likely to mimic these signatures. Also, at these energies 
there is negligible intrinsic v r background in the con- 
ventional atmospheric neutrino flux [15]. The prompt 
v T flux from charm hadron decays in cosmic-ray- induced 
air showers is also expected to be small [16Hlq . The 
signal v r are expected to come from the vicinity of the 
horizon since there is insufficient material for interactions 
in the downward-going direction and v T passing through 
the Earth emerge [19| at energies too low to create a UHE 
signature. 

The v T event topology depends on how much of the 
event is contained in the detector, the v T energy, and 
the composition of the r decay products. In this work 
only non-muonic r decays were considered. A partially 
contained UHE v T having only the decay vertex of r in 
the fiducial volume is denoted a "lollipop," while one 
having only the production vertex of the r in the fidu- 



cial volume is denoted an "inverted lollipop." A fully 
contained UHE v r having both production and decay 
vertices well separated in the fiducial volume is denoted 
a "double bang" [2(| . Fig. [T] shows a simulated double 
bang event in the 22-string configuration of the IceCube 
detector (IC22). 

Applying criteria to identify lollipop, inverted lollipop 
and double bang signatures produced by v T interactions, 
we derived limits on the diffuse UHE v T neutrino flux. 
We assumed a flux ratio of u^.v^Mt = 1:1:1 for this anal- 
ysis. We used 282.4 live-days of data collected in 2007- 
2008 by IC22. We describe the IC22 detector in Sec- 
tion [II] and the experimental and simulated data samples 
in Section Hill We present our analysis in Section HVl and 
the results in Section [Vl We discuss systematic errors in 
Section IVT1 and our conclusions in Section IVFll 



II. THE ICECUBE 22-STRING DETECTOR 

The 22-string configuration of IceCube (IC22) was de- 
ployed in early 2007, began taking physics-quality data 
in May of that year, and ended at the transition to Ice- 
Cube's 40-string configuration in April 2008. Each string 
consists of 60 digital optical modules (DOMs) buried 
deep in the icecap at the South Pole, with regular 17 m 
vertical spacing from 1450 to 2450 m below the surface, 
for a total of 1320 DOMs. The strings are situated on 
a regular grid with 125 m horizontal interstring spacing, 
covering the area shown in Fig. [2j Each DOM houses 

IC22 Geometry 
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FIG. 2. (Color online) Top view of the IceCube 22 string de- 
tector. Each string is represented as a dot with its associated 
string number. 

a photomultiplier tube (PMT) to detect the Cherenkov 
light, electronics for pulse digitization and other func- 
tions, and remotely-controllable calibration light sources. 
To reduce the impact of PMT signals due to random 
noise, only detected signals with minimum 0.25 single 
photoelectron (pe) PMT pulse height were digitized by 
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two types of waveform digitizers in situ: the ATWD 
(Analog Transient Waveform Digitizer) and an fADC 
(Fast Analog to Digital Converter). The time resolu- 
tion of the ATWD (fADC) is about 3.33 ns (25 ns) with 
a readout time window of about 450 ns (6.4 /zs). Thus 
the ATWD is used to capture detailed waveform infor- 
mation on a short time scale while the fADC records less 
detailed information on a longer time scale. The ATWD 
also supports three channels with different gains (xl6, 
x2, and x0.25) to extend its effective dynamic range. 

To further remove random noise the digitized signal 
in a DOM was required to be in close temporal coin- 
cidence with a signal in neighboring DOMs. The sig- 
nals satisfying such a temporal condition in hardware are 
called LC (Local Coincidence) hits. LC hits were then 
checked to see whether or not they satisfied a software- 
based trigger that selected for patterns potentially pro- 
duced by a charged particle. Groups of hits that satisfied 
a trigger condition were packaged into "events." Higher- 
level "filter" algorithms were applied to each event, and 
those events passing one or more filter conditions were 
transmitted over satellite to the northern hemisphere for 
higher-level analysis. However, all the data satisfying 
the software trigger conditions were stored on tape and 
shipped to the northern hemisphere. The software trigger 
and filter conditions applied to the data used in this anal- 
ysis are described in the section below. For more detail 
on the design, construction and performance of IceCube 
in general, see [2~lTj24T ] and references therein. 



III. DATA 
A. Experimental Data 

The DOM signals satisfying the LC condition were re- 
quired by the online data acquisition (DAQ) system at 
the surface computing system in the IceCube Laboratory 
to satisfy a "simple majority trigger" condition under 
which eight or more DOMs reported signals in a 5 /is time 
window ("SMT8"). The IC22 trigger rate of 500 to 620 
Hz followed the seasonal variation in the cosmic-ray muon 
flux. The DAQ system grouped together DOM hits sat- 
isfying the trigger condition into an event using a broad- 
ened ±8 /is time window. Triggered events used in this 
analysis were accepted if they also satisfied the extremely 
high energy (EHE) filter applied to the data online at the 
South Pole to reduce low energy events consistent with 
background. The EHE filter required > 80 DOMs regis- 
tering hits in the event. 

We split off about 30% of the full IC22 dataset (82.4 
live-days, uniformly distributed in time across the data- 
taking period) to use in conjunction with simulated data 
in the design of our subsequent selection criteria. In keep- 
ing with our procedures for maintaining blindness in the 
analysis of data, and thereby reducing human bias in 
the analysis of the data, the final result is based on the 
application of these selection criteria, unaltered, to the 



remaining 70% of the dataset (200 live-days). 



B. Simulated Data 



We employed simulated data to develop criteria that 
enhanced a possible astrophysical neutrino signal while 
diminishing backgrounds from atmospheric neutrinos 
and cosmic-ray muons. Exclusive use of simulated data 
also permitted us to maintain blindness. For the signal, 
the ANIS (All Neutrino Interaction Simulation) pack- 
age [25| was used to produce each neutrino flavor sep- 
arately. They were generated following an £ _1 energy 
spectrum to enhance event statistics at higher energy 
where this analysis is sensitive. The events were then 
run through the IceCube detector simulation. The muon 
(electron) neutrinos were generated over all zenith angles 
in the energy range between 10(50) GeV to 10 EeV while 
tau neutrinos were generated between 1 TeV and 1 EeV. 

Cosmic-ray muon backgrounds were simulated by gen- 
erating air shower events using the CORSIKA pack- 
age [26j , then propagating the muons to and through the 
detector volume with the MMC package [27[ , and finally 
applying the detector simulation to the resulting set of 
particles. 

For solitary air showers, a two-component model (28j 
was used. In this model, the entire mass spectrum of cos- 
mic rays is approximated by only proton and iron compo- 
nents. Compared to Horandel's polygonato model |29| . 
the two-component model agrees better with experimen- 
tal data at higher energy (beyond 100 TeV) where this 
analysis is sensitive. The cosmic ray primaries are sam- 
pled with an E -2 spectrum. In this way we were able 
to produce events more efficiently at the higher primary 
energies that contribute most strongly to the background 
at ultra-high energies. The cosmic ray flux was then re- 
weighted to match the expected spectrum. 

The acceptance of IC22 admits the possibility of 
detecting muons from multiple quasi-simultaneous air 
shower events, so we also simulated muons from two 
coincident air shower events. (Higher multiplicities oc- 
cur at a negligible rate in IC22 and were not simu- 
lated.) For coincident air showers, Horandel's polygo- 
nato model of cosmic rays was used. Solitary(coincident) 
atmospheric air showers were generated with energies be- 
tween 10(0.6) TeV-100 EeV and zenith angles between 
0-90°. 

After event generation and detector simulation, the 
simulated data were processed in the same way as real 
data, i.e., with sequential applications of trigger and filter 
conditions, as described earlier. 
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IV. TAU NEUTRINO IDENTIFICATION 

A. Selection Criteria 

Based on the characteristics of simulated data, we for- 
mulated several event selection criteria to exploit the 
UHE v T signatures of a track plus one or two showers, in 
contrast to conventional pure track-like or pure shower- 
like events. Two such criteria use the reconstructed num- 
ber of photoelectrons (N pc ) per DOM as a function of 
time, denoted ND° M (t). Fig. [3] shows ND° M (t) for a 
simulated inverted lollipop (top) and a simulated muon 
event (bottom). Note that the times of the hits are with 
respect to the event trigger time which has an extended 
readout time window of ± 8 fis in IC22. For this reason, 
all the hit times exhibit at least a ± 8 (is offset. 
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FIG. 3. The quantity Np° M (t) for a simulated inverted lol- 
lipop (top) and a simulated muon event (bottom), with pri- 
mary particle energies of 25.4 PeV and 2.38 PeV, respectively. 

To exploit the power of ND° M (t), we devised a pa- 
rameter called "maximum current ratio" (IR max ), de- 
fined as the maximum of Ii n /I ut where I\, 



AQ i: 



(out 



/AT], 



(out)- 



n(out) 

Here, AQ in was the charge, mea- 



sured in photoelectrons, collected by the DOMs in a slid- 
ing time window of length ATj n . The time window was 



optimized in this analysis to be 1.2 /xs long. The corre- 
sponding "out" variables were the charge and time mea- 
sured outside the sliding time window (see Fig. |4|). As 
shown in Fig. [5J IR ma x is small for track-like events and 
large for events containing showers, such as those pro- 
duced by v r . Since the IR max cut is related to energy, 
it will be applied to data as the last cut together with 
the other energy related cut explained at the end of this 
sub-section. 
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FIG. 4. (Color online) Illustration of the IR max parameter. 
See text for the explanation. 




FIG. 5. (Color online) The parameter IRmax for simulated 
signal (green histogram for lollipop and blue points for all 
v T events) and background (red points for atmospheric muon) 
events, and for data (gray histogram) passing the EHE fil- 
ter. The separation of signal from background is evident. 
The IRmax distributions of inverted-lollipop and double bang 
events are also well-separated from the background. 

Although IRmax is very effective at distinguishing 
most simple track-like background events from signal 
events, highly energetic muons can stochastically deposit 
large amounts of energy along their track lengths via 
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Bremsstrahlung, pair production, or photonuclear inter- 
actions, potentially mimicking v T events. Fig.[5]shows an 
example of simulated muon with such a Bremsstrahlung 
whose IR max value could be similar to that of a i/ T . 
Theoretically, v T events are most likely to have a large 
Np OM (t) at one or both of the temporal edges of the 
event. In practice, v T events had a large Np ( ? M (t) in 
the earliest third due to the presence of highly scattered 
photons that extended the temporal edge of the event 
to much later times. We expect future analyses to be 
able to devise criteria that reduce the impact of these 
scattered photons. The "local charge density" (p q ) pa- 
rameter was introduced to remove events consistent with 
a large energy deposit away from either temporal edge. 
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FIG. 6. The quantity Npv M (t) for a simulated muon event 
(336 PeV) with an energetic Bremsstrahlung. 

Dividing each event into equal thirds in time, p q was 
defined in each third as the total charge-to-time ratio. 
The event was rejected when p q < 5 N po /ns in either the 
first or last third. Figure [7] illustrates the p q parameter 
and shows how it can distinguish v r events from muon 
Bremsstrahlung events. Figure |S] quantifies the ability of 
p q to separate signal from background. 

Additional selection criteria were applied to further re- 
move backgrounds. The flux of downward-going muons 
from cosmic-ray air shower events was reduced by imple- 
menting a "veto layer" in software, removing any events 
in which the average Z position of the first 4 hits (^i n it), 
was in the top 50 m of the detection volume. Downward- 
going muons were further removed using the approximate 
event velocity Vz (m/ns), constructed from the differ- 
ence between the positions Z cog and Z- ln i t , divided by the 
difference in their respective times, i.e., T cog and T ul i t , 
where Z cog (T cog ) were the Z position (time) of the cen- 
ter of gravity of all hit DOMs. The times here are calcu- 
lated using the average time of the hits used to calculate 
the Z positions. We removed events consistent with a 
downward direction by requiring Vz < —0.1 m/ns. 

Background events arising from muon stochastic pro- 
cesses at or near the bottom of the detector, events whose 



FIG. 7. The quantity N? c OM (t) for a simulated lollipop (top) 
and double bang (bottom) event. Np<? M (t) is shown for 
a simulated inverted lollipop in Fig. [3] (top) for an atmo- 
spheric muon event in Fig. [6] The atmospheric muon, with a 
Bremsstrahlung roughly in the middle of its contained track 
length, would be rejected by the cut on p q described in the 
text, whereas the lollipop, inverted lollipop and double bang 
would not because the bulk of the detected light occurs suffi- 
ciently early in the event. 



muon tracks may go undetected, are removed by re- 
stricting our sample to events that were reasonably well- 
contained in the detector fiducial volume. We required 
the average depth position of all DOMs with signals to 
satisfy Z cog > —330 m (as measured from the center of 
the detector). 

We also applied a generic topological selection by cal- 
culating the eigenvalues of the tensor of inertia (Tol) 
of pulse amplitudes (instead of conventional mass) [30j 
from hit DOMs of each event and keeping only those 
events that tended towards sphericity. Perfectly spher- 
ical events will have three equal Tol eigenvalues, while 
perfectly track-like events will have one eigenvalue equal 
to zero. We therefore required that the ratio of smallest 
eigenvalue to the sum of all three eigenvalues was > 0.1. 

Remaining lower energy events were further reduced 
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FIG. 8. (Color online) The logarithm of the local charge den- 
sity parameter (p q ) for signal (green histogram for lollipop and 
blue points for all v T events) and background (red histogram 
for atmospheric muon) events, and for data (gray histogram) 
passing the EHE filter. The separation of signal from back- 
ground is evident. The p q distributions of inverted- lollipop 
and double bang events are also well-separated from the back- 
ground. 



in number by requiring a minimum N pe for each event. 
We required both IR max > 200 and log 10 N po > 4.2, the 
values of which were based on an optimization that is 
described in the following section. Figure [5] shows the 
distributions of these two selection criteria for simulated 
signal, simulated background, and 30% of the data, prior 
to the overall optimization of all the selection criteria. 



B. Optimization of Selection Criteria 

The final values for IR max and N po were optimized by 
minimizing the Model Rejection Factor (MRF) [3l| be- 
fore applying them to the full dataset. We varied the 
values of IR max and log 10 N pc as shown in Fig. [TUl find- 
ing a shallow minimum at MRF- 0.89. At this MRF, 
the expected all- flavor signal and background were 3.52 
and 0.81 events, respectively, using the Waxman-Bahcall 
upper bound [J| for the signal neutrino flux normaliza- 
tion with E -2 spectrum. The corresponding optimized 
values are IR max > 300 and log 10 N po > 4.0. However, in 
order to be conservative in the face of limited simulated 
event statistics, we chose instead to use IR max > 200 and 
log 10 N pc > 4.2, resulting in an MRF= 0.93 and expected 
all flavor signal and background event counts of 3.18 and 
0.60, respectively. 



erence purposes and are summarized in Table U for simu- 
lated signals and Table [XT] for simulated background and 
30% of the data sample. It is evident from Table U that 
this analysis, though designed to be sensitive primarily 
to UHE v T signals, also had appreciable sensitivity to 
UHE u c and signals. The final limit described below 
will therefore be applicable to all neutrino flavors. Fig- 
ures [TT1fT4l show the distribution of event rates (Hz) for 
each cut parameter for simulated signal and background, 
and a sample of IC22 data after application of the EHE 
filter (Fig. QT]) and sets of selection criteria SI (Fig. H"2"j). 
S2 (Fig.QS]), and S3 (Fig.E]). 

The event selection efficiency of the different cuts can 
be obtained from Fig. [TS] (top), where the simulated 
background and 30% of the data sample match well at 
each cut level. Figure [15] (bottom) shows the effective 
area A c g for each neutrino flavor after application of 
the SMT8 trigger condition and the full suite of selec- 
tion criteria. Using simulated signal, A G ff is defined by 
$„A e ffT = Ndct, where $„ is the neutrino flux prior to 
any propagation or interaction effects in the Earth, T is a 
length of time, and Ndct is the number of detected events. 

Figure [TC] shows the distribution of the true zenith 
angle from the simulation as a function of the primary 
v T energy for the events that passed all the selection cri- 
teria. As expected, most v T were from near the horizon, 
with the angular acceptance peaking at about 100° from 
vertical. 



V. RESULTS 

After unblinding the remaining 200 live-days of data 
and applying all the selection criteria, three events re- 
mained in the data sample. The predicted background 
from all simulated sources was 0.60 ± 0.19 events. The 
remaining data events are shown in Fig. 1171 

From a detailed study of these events, we deter- 
mined that one was consistent with light produced by 
an AMANDA optical module observed to emit light in- 
termittently (Fig. [T71 top). A second event was quali- 
tatively consistent with background from a nearly hori- 
zontal muon interacting near the bottom of the detector 
(Fig. [T71 middle). The third event had the character- 
istics of a neutrino- induced shower (Fig. [T71 bottom), 
and was also in the final sample of an independent IC22 
analysis that searched for shower- like signals [33j |. How- 
ever, we can not rule out this event as being produced 
by a cosmic-ray muon accompanied by a stochastic high- 
energy Bremsstrahlung. We have conservatively included 
all three events in the final sample in the derivation of 
the final result. 



C. Selection efficiency 



VI. SYSTEMATIC AND STATISTICAL ERRORS 



The event rates for the selection criteria described in 
Section llV Al were grouped into sets (EHE, Sl-4) for ref- 



The systematic and statistical errors in this analysis 
were obtained using signal and background simulations 
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TABLE I. Predicted signal event rates with statistical error after application of each set of selection criteria. LP, ILP and DB 
represent lollipop, inverted-lollipop and double-bang, respectively. For signal rates, the flux was normalized to the WB bound. 
(The first column provides labels for reference purposes in subsequent tables and figures.) 











MC simulation 






Set No. 


Selection Criteria 






Signal v (E~ 2 ) 










LP 


ILP 


DB 




v„ 








x 10" 9 [Hz] 


x 10" 9 [Hz] 


x 10" 9 [Hz] 


x 10" 8 [Hz] 


x 10" 8 [Hz] 


x 10" 8 [Hz] 


EHE 


NDOM > 80 


3.48 ± 0.11 


3.54 ± 0.09 


4.45 ± 0.16 


50.5 ± 0.5 


119 ± 2.2 


39.9 ± 0.7 


Sf 


p q > 5 N pc /ns 


3.42 ± 0.11 


3.05 ± 0.08 


4.30 ± 0.16 


24.0 ± 0.3 


29.3 ± 0.8 


23.9± 0.6 


S2 


Zait < 450 m, Vz < -0.1 m/ns 


2.55 ± 0.10 


2.91 ± 0.08 


3.95 ± 0.16 


22.6 ± 0.3 


24.9 ± 0.8 


22.9 ± 0.5 


S3 


Z cog > 330 m, Tol > 0.1 


2.32 ± 0.10 


2.29 ± 0.08 


3.02 ± 0.14 


15.7 ± 0.3 


11.8 ± 0.6 


17.5 ± 0.5 


S4 


IR max > 200, log 10 N pc > 4.2 


1.72 ± 0.08 


1.72 ± 0.06 


2.07 ± 0.11 


5.63 ± 0.08 


3.70 ± 0.15 


9.08 ± 0.2 



TABLE II. Predicted background event rates with statistical error after application of each set of selection criteria. For 
conventional neutrinos (labeled "conv" in the table), the Bartol model [l5l ] was used. For prompt neutrinos, the Martin GBW 
model [32] was used for v T and the Sarcevic standard model [TtJ was used for v^, and i/ e . 



Set No. 






MC simulation 






Data 






Background v 






Background /i 


30 % sample 




conv 


conv 


prompt 


prompt 
M 


^prompt 








x 10" 8 [Hz] 


x 10" 8 [Hz] 


x HT 10 [Hz] 


x 10" 8 [Hz] 


X 10" 8 [Hz] 


x 10~ 6 [Hz] 


x 10~ 6 [Hz] 


EHE 


184 ± 14.0 


6.88 ± 0.26 


33.4 ± 0.4 


23.6 ± 0.50 


9.95 ± 0.13 


830,000 


1,370,000 ± 438 


SI 


8.21 ± 1.80 


0.96 ± 0.06 


9.74 ± 0.17 


2.19 ± 0.12 


3.46 ± 0.05 


303 


246 ± 5.9 


S2 


8.11 ± 1.80 


0.96 ± 0.06 


9.62 ± 0.17 


2.05 ± 0.12 


3.42 ± 0.05 


41.2 


53.3 ± 2.7 


S3 


4.16 ± 0.66 


0.70 ± 0.06 


7.12 ± 0.14 


1.26 ± 0.09 


2.55 ± 0.04 


14.4 


20.8 ± 1.7 


S4 


0.24 ± 0.06 


0.04 ± 0.003 


0.91 ± 0.03 


0.15 ± 0.02 


0.43 ± 0.01 


0.026 ± 0.01 






and are summarized in Table IIIII In the following sub- 
sections, systematic errors on signal and background are 
explained followed by our result including both errors. 

A. Systematic Errors for Signal 

The systematic error due to our lack of precise knowl- 
edge of the DOM sensitivity to photons was obtained by 
simulating the effect of setting it to 90% and 110% of its 
nominal value resulting in [-4.7%, +7.9%] error. The sys- 
tematic error in the event rates reflecting uncertainties on 
the optical properties of the ice was obtained by simulat- 
ing events using different ice models. The ice models were 
created from data generated using in situ light sources. 
The baseline ice model 34] for this analysis used opti- 
cal properties of the ice measured at AMANDA depths 
and extrapolated to IceCube depths, while an alterna- 
tive ice model [35[ obtained them with a direct fit to the 
full range of IceCube light source data. Comparing the 
predictions of the two ice models resulted in a +29.4% 
error. 

The systematic uncertainty in the neutrino cross sec- 
tion came from two sources. One was from theoretical 
uncertainty in the parton distribution function evalua- 
tion and structure function and the other was from errors 
in the experimental measurement of the parton distribu- 
tion function by HERA [361 ] . From these two sources 



we estimated the systematic error in the neutrino cross 
section as ±6.4%. Very high energy events could satu- 
rate PMTs by exceeding the PMT's dynamic range. This 
could result in an incorrect estimation of the original neu- 
trino energy. Since the observable quantity most closely 
related to the energy is N pe , the systematic error associ- 
ated with the PMT saturation was obtained by observing 
the impact of changing the N pe cut from 90% to 110% 
of its original value. This error was found to be [-5.7%, 
+5.0%]. 

B. Systematic Errors for Background 

The systematic errors due to uncertainties in DOM 
sensitivity, ice properties, and DOM saturation behav- 
ior were obtained in the same manner as for the signals, 
as described in Section IVI AI They were estimated as 
[-4.7%, +7.9%], [-62%, +85%], and [-28.9%, +5.3%], re- 
spectively. 

In addition, there are systematic errors which applied 
only to the background. The muon event rate is known 
to change as a function of the atmospheric temperature 
above the South Pole plateau. Since our muon simula- 
tion assumed a rate pegged to that seen in October, the 
seasonal variation was taken into account as a systematic 
error and was estimated as [-24%, +18%] when compared 
with IC22 data at EHE filter level. The systematic er- 
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FIG. 9. (Color online) Distributions of the quantities IRmax vs. log 10 N pe for 30% of the data (row 1, left), simulated cosmic- 
ray background (row 1, right), simulated atmospheric neutrinos (row 2, left), simulated prompt atmospheric neutrinos (row 2, 
middle), simulated all-flavor neutrino signal (row 2, right), and v T inverted lollipop (row 3, left), lollipop (row 3, middle) and 
double bang (row 3, right), assuming an E~ 2 spectrum and prior to selection criteria optimization. The color code represents 
the event rate in Hz except for the data where it represents number of events in 30% of the data sample (82.4 live-days). The 
region in the upper right of each plot, indicated by the black lines, designates the region selected by the optimized criteria as 
described in Section [IV Bl 
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• IRmax >= 


50 


• IR m ax >= 


100 


• IRpnax >= 


150 


• IRmax >= 


200 


• IRmax >= 
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FIG. 10. (Color online) We tuned the selection criteria 
IRmax and log 10 N pe using the Model Rejection Factor (MRF) 
formalism. The plot shows how the MRF varies as a function 
of log 10 N pc for different values of IRmax- We chose values for 
these parameters near but not exactly at the minimum shown 
for reasons explained in the text. 
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FIG. 12. (Color online) Distributions of selection criteria for 
signal, BG, and 30% of the IC22 data just before application 
of the "S2" set of selection criteria. The vertical (horizontal) 
arrow lines represent the cut values (selected regions) of that 
set. 



FIG. 11. (Color online) Distributions of selection criteria for 
signal, BG, and 30% of the IC22 data just before application 
of the "SI" set of selection criteria. The vertical (horizontal) 
arrow line represents the cut value (selected region) of that 
set. 



ror due to cosmic ray composition was also obtained by 
switching constants and slopes between proton and iron 



in the two component model data. At S3, just before the 
final cut to have enough statistics, we obtained -24% by 
this method. 

There are alternative models for the prompt neutrino 
flux. For this analysis, the base models used for the 
prompt neutrino flux are Sarcevic standard flux model 
for Vu and v e [l7[, and Martin GBW model for v T [32l |. 
As an alternative, we have also considered the Sarcevic 
minimum and maximum flux models 17], from which we 
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FIG. 13. (Color online) Distributions of selection criteria for 
signal, BG, and 30% of the IC22 data just before application 
of the "S3" set of selection criteria. The vertical (horizontal) 
arrow lines represent the cut values (selected regions) of that 
set. 



estimate a [-59%, +30%] systematic error on the prompt 
neutrino flux. 



C. Result including Statistical and Systematic 
Errors 



Since it was computationally feasible to generate a 
large amount of simulated signal, the statistical error 
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FIG. 14. (Color online) Distributions of selection criteria for 
signal, BG, and 30% of the IC22 data just before application 
of the "S4" set of selection criteria. The vertical (horizontal) 
arrow lines represent the cut values (selected regions) of that 
set. The legend for these plots is the same as in Fig. [13] 



on the simulated signal is small (±2.3%). By contrast, 
the considerably larger statistical error on the simulated 
background (±39%) reflects the aggregate effect of the 
high rejection efficiency of our selection criteria and the 
limitations imposed by finite computational resources. 

We combined the systematic and statistical errors 
and applied a profile log-likelihood method [37| to de- 
rive the following limit on the astrophysical all-flavor 



90% CL 

-1 o-l 



E 2 <f>{v K ) < 16.2 



neutrino flux at 

10 -8 GeV cm~ 2 sr _i s _i for the 3 observed events from 
the 200 live-days of IC22 data. This limit applies to 
the primary neutrino energy range of 340 TeV < E v < 
200 PeV, covering the middle 90% of the accepted sim- 
ulated signal. Fig. [TH] shows this limits together with 
several theoretical model predictions. The upper limit 
on the tau neutrino flux is one third that of the all-flavor 
astrophysical neutrino flux if one assumes a flavor ratio 
of Vc.v^-.Vr = 1:1:1. 
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FIG. 15. (Color online) Top: Event rate (Hz) at each cut level 
for simulated signal and background, and 30% of the data 
sample. At S4, there were zero events in the data sample, 
and so the 90% CL upper limit value was plotted as indicated 
by the black arrow. Bottom: IC22 effective areas for each 
neutrino flavor (black: v^, red: v c , blue: u T ) after application 
of the SMT8 trigger (dashed lines) and after application of all 
selection criteria (solid lines). 



VII. CONCLUSIONS AND OUTLOOK 

A set of selection criteria designed for UHE 
u T detection were applied to IceCube data. These cri- 
teria also had appreciable efficiency for UHE v c and 
ffj, detection. We applied these criteria to 200 live-days of 
data from IceCube's 22-string configuration and observed 
3 events in the final sample. We therefore set a 90% 
CL upper limit on the astrophysical UHE all-flavor neu- 
trino flux of E 2 $(v x ) < 16.2 x 10" 8 GeVcm" 2 sr" 1 s" 1 . 
The analysis improves on the previous limit set by 
AMANDA [H [li, SH with comparable integrated ex- 
posure. Future IceCube searches specialized for v r will 
be more sensitive due to the increased fiducial volume 
relative to IC22. The large volume will also warrant the 
application of sophisticated v r reconstructions, further 
improving the sensitivity of these searches. 




FIG. 16. (Color online) The distribution of the true zenith an- 
gle from the simulation as a function of the primary v r energy 
for the events passing all selection criteria. Color code repre- 
sents event rate in Hz. 
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TABLE III. Summary of the systematic and statistical errors 
for signal and background events from the simulated data. 



Source 



DOM sensitivity 
Ice properties 
v cross section 
PMT saturation 
Cosmic ray flux 
Cosmic ray composition 
Seasonal variation 
Prompt v flux model 



Total syst. error 



Signal 



4.7%, + 7.9% 
- 0%, + 29% 
6.4%, + 6.4% 
5.7%, + 5.0% 

N/A 

N/A 

N/A 

N/A 



7.9%, + 31% 



Background 



4.7%, + 7.9% 

- 62%, + 85% 

N/A 

- 29%, + 5.3% 

- 0%, + 16% 

- 24%, + 0% 

- 24%, + 18% 

- 59%, + 30% 



97%, + 94% 



* 





















] 


1 


1 









FIG. 17. (Color online) Diagrams of the three events surviv- 
ing the final selection criteria applied to the 200 live-days of 
IC22 data. The radius of each circle is proportional to the 
number of photons detected by the PMT. 
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FIG. 18. (Color online) The limits on production of UHE 
neutrinos. JA]: AMANDA-II cascade all flavor limit (1001 
live-days) [3E|, [B]: IC22 cascade all flavor limit (257 live- 
days) [33|], [C]: Baikal all flavor limit (1038 live-days) 4011 . 

[D] : AMANDA-II UHE all flavor limit (457 live-days) [4l(, 

[E] : IC22 UHE all flavor limit from diffuse sources using the 
analysis described in this paper (this work, 200 live-days), 
[E*]: IC22 UHE all flavor sensitivity (this work, 200 live- 
days), [F]: ANTARES '07-'09 v„ x 3 334 d [H, [G]: IC40 

x 3 (375.5 live-days) [aa]: Waxman Bahcall (v^ and 
TV ) model 1998 x| Q|, [bb]: Stecker AGN (Seyfert) 2005 0], 
[cc]: Waxman Bahcall Prompt GRB model , and [dd]: At- 
mospheric neutrino flux (Bartol + Sarcevic standard model) 
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